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We study the effect of two-dimensionality on step bunching induced by drift of adatoms. When 
anisotropy of the diffusion coefficient changes alternately on consecutive terraces like a Si(001) vicinal 
face, bunching occurs with the drift of adatoms. If the fluctuation of step bunches is neglected as 
in the one- dimensional model, the bunching with step-down drift is faster than that with step-up 
drift in contradiction with the experiment by Latyshev and coworkers. In a two-dimensional model, 
the step bunches wander heavily with step-up drift and recombination with neighboring bunches 
occur more frequently than those with step-down drift and the bunching is accelerated. When the 
difference of kinetic coefficients between two types of steps is taken into account, the bunching with 
step- up drift can be faster than that with step-down drift. 
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I. INTRODUCTION 

The Si(001) surface is reconstructed by dimcrization of 
surface atoms. When its vicinal face is tilted in the (110) 
direction, the terraces with dimer rows parallel to the 
steps (we call Ta) and those with dimer rows perpendic- 
ular to the steps (we call Tb) appear alternately. Since 
the surface diffusion along the dimer rows is faster than 
that perpendicular to the dimer rows, the anisotropy of 
the surface diffusion changes alternately on consecutive 
terraces. 

On the vicinal face, two types of step instabilities, step 
wandering and step bunching, occur when a specimen is 
heated by direct electric current Q, • The step wander- 
ing occurs with step-up current in the region of relatively 
large inclination (the tilting angle is 0.08° < 6 < 0.5°) @. 
Due to the step wandering, grooves perpendicular to the 
steps appear on the vicinal face. The step bunching oc- 
curs irrespective of the current direction in the region of 
small inclination (6 < 0.08°). The type of dominant ter- 
races, which separate step bunches, is Tb with step-down 
current and Ta with step-up current. The size of bunch 
increases with time as t 1 / 2 0, which is independent of 
the drift direction. The growth rate of the bunches with 
step-down current seems slightly slower than that with 
step- up current. 

The step instabilities are caused by drift of adatoms 
induced by the current. By taking account of alternation 
of the anisotropic surface diffusion, the step instabilities 
are theoretically explained. If the repulsive interaction 
is strong so that the step bunching is suppressed, the 
step wandering occurs with step-up drift £| , and straight 
grooves parallel to the drift appears by in-phase wander- 
ing. The motion of the steps is given by the solution of 
the nonlinear equation derived by Pierre-Louis and co- 
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workers [f| . The step bunchingoccurs irrespective of the 
direction of the drift H 0, HH, H E3 Since the cur- 
rent and the drift are in the same direction II4 llfl , the 
results are consistent with the experiments jjj, |3, LJ ■ In a 
one-dimensional step flow model [ljj , the size of bunches 
increases with time as t$ with j3 slightly smaller than 
1/2, consistent with the experiment Q. However, the 
growth rate with step-up drift is slower than that with 
step-down drift since terraces with fast diffusion in in- 
direction are dominant with step-down drift. This is in 
contradiction with the experiment Q- Since the model 
with alternating diffusion anisotropy has explained the 
bunching and the wandering instabilities on the Si(001) 
vicinal face consistently, this disagreement is a major ob- 
stacle to the unified understanding. 

In the one-dimensional step flow models [ToL ITll ] , the 
motion of step bunches with step-down drift is similar 
to that with step-up drift except the time scale. In 
the Monte Carlo simulation, however, the step pattern 
is changed by the drift direction 0, Q : the step bunches 
with step-up drift wander more than that with step- 
down drift. Such a difference of the step motion in 
two-dimension may solve the disagreement in the growth 
rate between the experiment and the one-dimensional 
model |lfjj . In this paper, we carry out Monte Carlo sim- 
ulations and show that the difference of the growth rate 
vanishes in the two-dimensional model. 



II. MODEL 

For simplicity, we use a square lattice model with the 
lattice constant a = 1. We take x-axis parallel to the 
steps and y-axis in the down-hill direction. Boundary 
conditions are periodic in the x-direction and helical in 
the y-direction. We forbid two-dimensional nucleation 
and use solid-on-solid steps, i.e. the step positions are 
single valued functions of x. 

We repetitively select a solid atom at the step or an 
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adatom on the terrace. We perform the diffusion and 
solidification trial for the adatom and melting trial for 
solid atom. In the diffusion trial, the adatom hops to a 
neighboring site. The anisotropy of the diffusion coeffi- 
cient and the drift of adatoms are taken into account in 
the hopping probability. On Ta, where the surface dif- 
fusion in the x-direction is faster, an adatom on the site 
moves to (i ± 1, j) with the probability 1/4 and to 
(i, j ± 1) with the probability Pd(l ± Fa/k B T)/4, where 
Pd(< 1) is the ratio of the two diffusion coefficients. F is 
the force to cause the drift. F > represents the drift 
in the down-hill direction. On Tb, where the surface dif- 
fusion in the y-direction is faster, an adatom on the site 
moves to (i ± l,j) with the probability pd/4 and 
to (i, j ± 1) with the probability (1 ± Fa/k B T)/4. For a 
diffusion trial, the time increment is At — l/4iV a , where 
N a is the number of adatoms so that the fast diffusion 
coefficient is unity. 

If the adatom comes in contact with a step from the 
lower terrace after a diffusion trial, solidification occurs 
with the probability 
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AE S is given by AE S = ex (the increment of the step 
perimeter) and <fi is decrement of the chemical potential 
by solidification. The step stiffness $ is related to e as 
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If we select a solid atom, melting trial is performed. 
When an adatom is absent on the top of the solid atom, 
melting occurs with the probability 



Pm = 



1 + exp 
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There is no extra diffusion barrier over the steps: steps 
are permeable. 



III. RESULTS OF SIMULATION 

In our simulations, we neglect the long-range repulsive 
interaction between steps, but take into account a short- 
range repulsive interaction by forbidding overlap of steps. 
Impingement of atoms and evaporation are absent. 

Figures n and [5| represent snapshots of the step bunch- 
ing. System size is 256 x 256 and the number of steps 
is 64. The parameters are e/k B T = 1.0, <j>/k B T = 1.5, 
Fa/k B T = ±0.08 and pd = 0.5. Initially a few adatoms 
are present on the vicinal face. The dotted lines represent 
Sa steps and the solid lines represent Sb steps. 

The vicinal face is unstable with the drift of adatoms. 
Pairing of Sa and Sb occurs in the initial stage. The 
upper side step in a pair is Sa with step-up drift and 
Sb with step-down drift. Small bunches arc formed by 
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FIG. 1: Snapshots of step bunching at t — 7.1 x 10 2 (a) with 
step-down drift and (b) with step-up drift 



coalescence of step pairs. Since the stiffness is small, the 
bunches wander and connect with each other at many 
places (Fig.^). From the figures one may have impression 
that he step with step-down drift are more straight. 
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FIG. 2: Snapshots of step bunching (a) with step-down drift 
at t = 3.6 x 10 4 and (b) with step-up drift at ( = 3.5x 10 4 . 

The effect of the drift direction on the form of bunches 
becomes evident in a late stage (Fig.EJ. With step-down 
drift, the bunches are straight and there are few recom- 
bination of bunches. With step-up drift, the wandering 
width of the bunches is large. The bunches collide with 
each other and frequent recombination is seen. The dif- 
ference of the form may affect the time evolution of bunch 
size. 

To test the effect of the wandering and recombina- 
tion on the growth rate, we carry out simulations with 
a narrow system (Fig. The number of steps is 128 
and the system size is 16 x 512. We use the parameter 
e/k B T = 0.5 for step bunches to wander easily. Other 
parameters are the same as that in Fig. [21 The wan- 
dering of step bunches with step-up drift is suppressed 
because of the narrow system width. The step bunches 
arc straight irrespective of the drift direction. 

In Figure |3| the number iV max of steps in the largest 
bunch at x = 1 is plotted as a function of time. The 
bunches grow by the collisions of straight bunches due 
to the fluctuation of position of bunches [lfl E3. The 
growth rate with step-up drift is slower than that with 
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with r = 0.1. Other parameters and system size are the 
same as those in Fig. The growth exponent does not 
change significantly. The growth rate with step-down 
drift is suppressed as expected from the reduction factor 
r, but that with step-up drift is slightly enhanced. 
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FIG. 3: Time evolution of the largest bunch size with step-up 
drift (A) and step-down drift (O). The system size is 16 x 512 
and the number of steps is 128. The result is obtained by 
averaging over 20 runs. 



step-down drift like the one-dimensional model. A^ max 
seems to increase in the power law, N max ~ t@ with 
(3 « 0.4 ^{|, which is consistent with the one-dimensional 
result with a combination of lnr and r~ 2 potentials [ill ]. 

We carry out simulations with a large system size: 
512 x 512 (Fig. 0}. In contrast to the one-dimensional 
model 1 1 1| . the growth of bunch size with step-up drift 
is as fast as that with step-down drift. Thus the slow 
diffusion in y-direction is compensated by the efficient 
coalescence by the wandering of bunches. 
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FIG. 4: Time evolution of the largest bunch size with step-up 
drift (A) and step-down drift (o). The system size is 512 x 512 
and the number of steps is 128. The result is obtained by 
averaging over 20 runs. 

In the experiment 0, however, the step bunching with 
step-up current seems slightly faster than that with step- 
down current. Additional mechanism to accelerate the 
bunching with step-up drift may be required. On the 
Si(001) vicinal face, Sb steps are rougher than Sa steps 
and the kinetic coefficient is probably larger for Sb- We 
take account of the difference as the probabilities of so- 
lidification and melting. At Sa, the probability of solidi- 
fication is assumed to be rp s and that of melting is rp m 
with r < 1. 

Figure [S] represents the time evolution of bunch size 
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FIG. 5: Time evolution of the largest bunch size with step-up 
drift (A) and step-down drift (o). The step kinetics of Sa is 
slower (r = 0.1) The system size is 512 x 512 and the number 
of steps is 128. The result is obtained by averaging over 30 
runs. 

In figure we show that the step patterns. As ex- 
plained already, the bunches with step-down drift are 
straight. The bunches with step- up drift are wavy and 
shows more recombination patterns. The reduction of 
kinetic coefficient of Sa seems to enhance the waviness. 
As a result, the growth rate with step-up drift becomes 
faster than that with step-down drift, which is in agree- 
ment with the experiment Q|. 
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FIG. 6: Snapshots of step bunching (a) with step-down drift 
at t = 3.6 x 10 4 and (b) with step-up drift at ( = 3.5x 10 4 . 



IV. SUMMARY AND DISCUSSION 

We studied the time evolution of bunch size by Monte 
Carlo simulation. With step-down drift, the step bunches 
are straight and there are few recombination of bunches. 
With step- up drift, the bunches wander and recombina- 
tions of bunches occur frequently 
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In the one-dimensional model [lOllllj . the velocity of a 
bunch is roughly proportional to the diffusion coefficient 
of y-direction in large terraces. Bunches with step-down 
drift grow faster than those with step-up drift, which is 
in disagreement with the experiment |2|. In the two- 
dimensional model, bunches with step-up drift wander 
and collide with each other more frequently than those 
with step-down drift. Since coalescence of bunches starts 
from the connected parts, the step bunching with step- 
up drift occurs more frequently than that with step-down 
drift. Also the fast transverse diffusion with step-up 
drift helps recombination of bunches. Consequently the 
growth rate of bunches with step-up drift is as fast as 
that with step-down drift. 

When the difference of the kinetic coefficient is taken 
into account, the bunches with step- up drift grows faster 



than those with step-down drift. The reversal of the 
growth rate is thus attributed to the two-dimensional 
step pattern. Since we have not succeeded in quanti- 
fying the effect, it is not clear if the proposed effect is 
sufficient to explain the experiment Q. 



Acknowledgments 

This work was supported by Grant-in- Aid for Scientific 
Research from Japan Society for the Promotion of Sci- 
ence. M. U. and Y. S. benefited from the inter-university 
cooperative research program of the Institute for Mate- 
rials Research, Tohoku University. 



[1] L. V. Litvin, A. B. Krasilnikov, and A. V. Latyshev, Surf. 

Sci. 244, L121 (1991). 
[2] A. V. Latyshev, L. V. Litvin, and A. L. Aseev, Surf. Sci. 

130-132, 139 (1998). 
[3] J.-F. Nielsen, M. S. Pettersen, and J. P. Pelz, Surf. Sci. 

480, 84 (2001). 
[4] M. Sato, M. Uwaha, Y. Saito and Y. Hirose, Phys. Rev. 

B 67, 125 (2003). 
[5] O. Pierre-Louis, C. Misbah, Y. Saito, J. Krug, and P. 

Politi, Phys. Rev. Lett. 80, 4221 (1998). 
[6] S. Stoyanov, Jpn. J. Appl. Phys. 29, L659 (1990). 
[7] A. Natori, H. Fujimura, and H. Yasunaga, Jpn. J. Appl. 

Phys. 92, 1164 (1992). 
[8] A. Natori, H. Fujimura, and M. Fukuda, Appl. Surf. Sci. 

60/61, 85 (1992). 
[9] M. Sato, M. Uwaha, and Y. Saito, J. Cryst. Growth 237- 

239, 43 (2002). 



[10] M. Sato, M. Uwaha, T. Mori, and Y. Hirose, J. Phys. 

Soc. Jpn 72, 2850 (2003). 
[11] M. Sato, T. Mori, M. Uwaha, and Y. Hirose, J. Phys. 

Soc. Jpn 73, 1827 (2004). 
[12] M. Ichikawa and T. Doi, Appl. Phys. Lett. 60, 1082 

(1992). 

[13] J.-J. Metois, J.-C. Heyraud, and A. Pimpinelli, Surf. Sci. 
420, 250 (1999). 

[14] M. Sato, M. Uwaha, Y. Saito, and Y. Hirose, Phys. Rev. 
B 67, 125408 (2003). 

[15] In the limited size of our simulation, it is difficult to define 
a bunch size. We regard a step pair as being contained in 
the bunch if the distance between the step pair and the 
bunch is smaller than a prescribed value. The obtained 
exponent, though not very sensitively, depends on the 
value. 



